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SUMMARY 

Consideration  is  given  to  the  bearing  ionospheric  scintill¬ 
ation  has  on  the  design  and  on  the  evaluation  of  the  perform¬ 
ance  of  satellite  communication  links.  Emphasis  is  placed  on 
establishing  a  design  parameter  of  universal  application. 

The  cumulative  amplitude  distribution  function  (cdf) ,  which 
has  found  application  in  this  regard,  is  shown  experimentally 
to  be  well  represented  by  the  Nakagami  m-distribution  at  all 
scintillation  levels.  However,  the  cdf  is  demonstrated  to  be 
strictly  valid  only  for  amplitude  modulated  signals  of  short 
duration.  Autocorrelation,  power  spectra  and  level  crossing 
techniques  are  investigated  in  an  effort  to  establish  a  design 
parameter  more  universally  applicable  than  the  cdf.  A 
quantity  message  reliability,  is  defined  which  removes  the 
short  duration  restriction  inherent  in  the  cdf.  A  formula 
for  message  reliability  is  developed  which  allows  it  to  be 
evaluated  for  any  fade  margin  and  message  length,  given  the 
scintillation  index  and  the  half  width  of  the  autocorrelogram 
of  the  scintillating  channel.  A  design  chart,  embodying  this 
information,  is  also  illustrated. 
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1 .  INTRODUCTION 

ConMunication  links  involving  transionospheric  propagation  are  subject  to  dis¬ 
ruption  due  to  ionospheric  scintillation(ref . 1) .  The  scintillation  takes  the 
form  of  random  fluctuations  in  both  the  amplitude  and  the  phase  of  the  received 
signal.  This  effect  is  most  marked  in  the  VHF/UHF  region  of  the  electromagnetic 
spectrum  but  it  is  also  evident  at  SHF  and  above (ref .2) .  The  effect  is  due  to 
perturbations  introduced  into  a  steady  incident  wave  by  irregularities  of  electron 
density  which  exist  mainly  in  the  F-layer  of  the  ionosphere (ref . 1) .  The 
occurrence  of  these  irregularities  has  been  studied  for  a  number  of  years (ref. 3) 
and  their  global  occurrence  morphology  has  been  documented  to  such  an  extent 
that  this  morphology  is  now  amenable  to  successful  model ling (ref .4,5,6) . 

However,  this  report  is  concerned  not  with  the  occurrence  morphology  but  rather 
with  the  nature  of  the  degradation  of  a  communications  channel  caused  by  the 
presence  of  scintillations. 

Scintillations  cause  both  enhancements  and  fading  of  the  signal  about  a  median 
level  as  well  as  random  fluctuations  in  phase.  When  scintillation  fades  occur 
which  exceed  the  fade  margin  built  into  a  communications  link,  the  performance 
of  the  link  will  be  degraded.  The  degree  of  degradation  will  depend  on  how  far 
the  signal  fades  relative  to  the  fade  margin,  the  duration  of  the  fade,  the  type 
of  modulation  employed  and  the  criteria  of  acceptability  for  the  channel.  In 
this  report  discussion  of  these  factors  will  be  illustrated  with  an  analysis  of 
UHF  data  obtained  with  the  aid  of  MARI SAT  II (ref. 7). 

The  report  commences  in  the  next  Section  by  considering  the  source  of  the  data 
employed.  This  is  followed  in  Section  3  by  theoretical  and  experimental 
considerations  of  the  effect  scintillations  have  on  signals  of  short  duration. 
Section  4  extends  these  considerations  to  signals  of  longer  duration.  Finally, 
the  conclusions  reached  are  summarized  and  discussed  in  Section  5. 


2.  THE  DATA 

The  scintillation  data  analysed  in  the  following  Sections  was  obtained  by 
observing  MARI SAT  II  at  Manus  Island  (147.37°E,  2. 04°S) (ref . 7) .  MARI SAT  II  is 
at  synchronous  altitude  over  the  equator  at  176. 5°E  and  hence  was  viewed  a  little 
north  of  east  and  at  an  elevation  of  about  56°  from  Manus  Island.  The  carrier 
at  256.550  MHz,  known  as  the  Fleet  Broadcast  Channel,  which  operates  with  time 
division  multiplexed  phase  shift  keyed  modulation  at  1200  bits/s  was  monitored 
with  a  receiver  of  high  dynamic  range.  An  IF  bandwidth  of  7  kHz  was  used  while  a 
post-detection  bandwidth  of  5  Hz  was  employed.  The  receiver's  AGC  voltage  was 
digitized  with  6  bit  sampling  and  the  digital  data  was  recorded  on  magnetic  tape. 
At  least  one  calibration  was  carried  out  each  day,  which  allowed  the  digitized 
AGC  data  to  be  converted  in  the  subsequent  processing  to  equivalent  signal  power 
(dBm)  at  the  antenna. 

In  digitizing  the  data  a  sample  rate  of  10  samples/s  was  used.  This 
corresponds  to  a  Nyquist  frequency  of  5  Hz.  Since  the  scintillation  frequencies 

of  interest  are  less  than  1  Hz,  the  resultant  oversampling  of  the  analog 
information  affords  some  protection  against  aliasing.  This  protection  was 
further  enhanced  by  the  action  of  the  post-detection  filter. 

It  is  planned  to  provide  the  receiving  equipment  with  the  capability  of 
recording  the  bit  error  rate  detected  in  the  phase  shift  keyed  modulation  of 
the  satellite's  transmissions.  The  bit  error  rate  is  expected  to  be  sensitive 
to  the  signal's  phase  stability  as  well  as  its  power  stability,  whereas  the  AGC 
measurements  are  a  measure  of  the  signal  power  stability  only.  However,  no 
bit  error  rate  data  were  available  at  the  time  this  report  was  written. 
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3.  SMALL  DURATION  SIGNALS 

In  this  Section  consideration  will  be  given  to  the  effect  scintillations  have 
on  communications  channels  which  are  carrying  signals  of  short  duration, 
e.g.  digital  data  at  rates  of  the  order  of  kilobits  per  second  and  higher. 
Consideration  of  the  effects  scintillations  have  on  signals  of  relatively  long 
duration,  e.g.  as  in  voice  communications,  will  be  deferred  to  Section  4. 

3 . 1  Theory 

In  order  to  quantify  the  fluctuations  in  the  amplitude  (R)  of  a 
scintillating  signal,  a  quantity  called  the  scintillation  index  (S4)  can  be 
defined  as  follows (ref .8) . 


s4  =  UR4  -  (R2)2] /(R2)2]%  (1) 

This  index  proves  to  be  a  convenient  measure  of  scintillation  activity 
when  the  observation  of  scintillations  is  used  as  a  technique  for  studying 
ionospheric  irregularities.  This  is  because  the  index  can  be  related 
directly  to  the  properties  of  the  irregularities(ref.8,9) .  However,  in  the 
engineering  of  communication  circuits  employing  transionospheric  propagation, 
a  more  system  oriented  measure  is  needed,  such  as  the  probability  of 
communication  being  lost  due  to  scintillation  effects. 

For  most  communication  systems,  the  "outage”  probability  can  be  defined 
in  terms  of  the  probability  of  the  signal  level  falling  below  some  designated 
value.  For  example,  in  the  interval  T  for  the  scintillating  signal 
illustrated  in  figure  1,  the  outage  probability  U(x0)  at  a  margin  of  Xo 
relative  to  the  median  level,  is  given  by 

T 

U(Xo)  =  ^  Stf/T 

t=0 


On  the  other  hand,  the  probability  of  the  signal  being  above  the  margin  Xo  is 


_T 

P(Xo)  =  ^  5ts/T  =  1  -  U(Xo)  (2) 

t=0 

The  signal  fluctuates  at  the  receiver  because  it  is  a  random  superposition 
of  random  vectorial  elements.  Nakagami (ref . 10)  has  investigated  such 
systems  extensively  and  has  shown  that  the  probability  distribution  p(R)  of 
the  resultant  signal  amplitude  (R)  is  represented  by 


P(R) 


„  m  r.2m-l 

2m  R_ 

r(m)  (R^)2 


This  is  known  as  the  m-distribution  and  it  can  be  demonstrated  that 


m 


i/sl 


(3) 
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It  is  worth  noting  that  the  m-distribution  reduces  to  the  Rayleigh 
distribution  in  the  special  case  where  m  =  1.  The  Rayleigh  distribution  is 
the  classical  result  for  the  distribution  expected  from  a  superposition  of 
random  vectors.  Strictly  it  applies  only  to  the  case  where  there  is  no 
steady  signal  component,  whereas  the  m-distribution  encompasses  both  this 
possibility  and  that  where  there  is  a  background  signal. 

Defining  the  signal  level  X  as 

X  =  10  logi  o  (R2  /R2 ) 

the  m-distribution  in  terms  of  signal  level  becomes 


P*  =  M  r(m)  exPt  2X/M  -  exp(2x/M)j] 


(4) 


where 


M  =  20  logic  e  =  8.686  (5) 

It  follows  that  the  probability  of  the  signal  level  falling  below  some 
specified  level  (i.e.  the  outage 

,  Xo 

U(Xo)  =  /  p  (X)dX 

j  _oo 


Similarly  the  probability  of  the 

CO 

P(Xo)  =  j  p'  (X)dx 
•>  Xo 

P (Xo )  has  become  known  as  the  cumulative  amplitude  probability  distribution 
function  or  cdf (ref. 11)  and  is  plotted  (as  a  percentage)  for  a  series  of 
values  of  m  in  figure  2.  Here  the  signal  levels  (Xo)  are  referred  to  the 
median  level  at  which  Xo  is  set  zero. 

For  amplitude  modulated  signals  of  sufficiently  short  duration,  the  out¬ 
age  rate  will  closely  approximate  the  instantaneous  outage  probability  U(Xo)- 
Hence  for  amplitude  modulated  signals  with  bit  rates  in  excess  of  about 
1  kbits/s,  the  bit  error  rate  (BER)  will  equal  U(x0).  That  is,  for  such 
signals 


probability)  is 


Xo  2mm 


M 


exp[m[2x/M-exp(2x/M)j]dX 


_oo 


signal  being  above  the  margin  Xo  is 


00  0  m 

2m 


=  J  M'f(m)  exp[m)2x/M-exp(2x/M)]]dX 


Xo 


(6) 


BER  =  100  -  cdf  (7) 

Bits  may  be  lost  from  signals  which  involve  phase  shift  keying  either 
because  of  amplitude  scintillation  or  because  of  phase  scintillation.  Hence 
for  such  phase  modulated  signals 


BER  >  100  -  cdf 


(8) 
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Because  of  the  greater  instrumental  difficulties  involved  in  their 
observation,  little  data  is  available  on  phase  scintillation.  However,  the 
random-phase  screen  theory  of  scintillation(ref . 12)  shows  that  the  spatial 
scale  of  the  amplitude  fluctuations  is  always  less  than  the  spatial  scale  of 


TABLE  1.  PARAMETERS  FOR  THE  PERIODS  OF  SCINTILLATION  ANALYSED 


Year 

Day 

Start 

L.T. 

s4 

m 

7  0. 5 
(s) 

f 

c 

(Hz) 

P 

1976 

237 

0030 

0.20 

24.9 

5.5 

0.025 

2 

to 

3 

1976 

236 

2330 

0.37 

7.2 

3.35 

0.04 

3 

to 

4 

1976 

236 

2300 

0.65 

2.3 

2.6 

0.05 

4 

to 

5 

1976 

236 

1930 

0.86 

1.3 

7.25 

(0.02) 

3 

to 

4 

1976 

236 

2000 

1.08 

0.9 

2.9 

(0.02) 

4 

to 

5 

the  associated  phase  fluctuations  in  the  diffraction  pattern  produced  by  a 
screen(ref . 13, 14) .  Further,  the  differences  in  the  scale  sizes  increases 

with  the  r.m.s.  phase  deviation  and  hence  with  scintillation  index.  Since 
the  amplitude  and  the  phase  scintillations  are  both  the  result  of  the  move¬ 
ment  of  the  diffraction  pattern  across  the  observing  point,  it  follows  that 
the  phase  fluctuations  are  generally  slower  than  the  associated  amplitude 
fluctuations,  the  disparity  increasing  with  increasing  scintillation  index. 

This  being  the  case,  amplitude  scintillations  will  be  the  predominant  factor 
in  determining  the  bit  error  rate  for  phase  modulated  signals  especially 
when  the  scintillation  index  is  large.  This  conclusion  is  apparently  borne 
out  by  observation (ref . 15) . 

3.2  Experimental  verification 

Because  of  the  connection  between  BER  and  cdf  noted  in  the  previous  Section, 
it  is  important  to  verify  that  the  m-distribution  is  a  valid  representation 
of  the  instantaneous  probability  distribution  of  the  signal  amplitude. 

While  verifications  of  this  type  have  been  carried  out  for  values  of  m  near 
unity (ref . 11, 16) ,  the  data  described  in  Section  2  allow  the  validity  of  the 
m-distribution  representation  to  be  evaluated  for  a  whole  range  of  values  of 
m. 

Figure  3(a)  shows  the  analog  representation  of  the  signal  strength 
fluctuations  during  a  period  of  30  min  commencing  on  0030  hours  L.T.  on  day 
237  of  1976.  During  this  period  the  application  of  equation  (1)  to  the 
digital  data  yields  a  scintillation  index  of  0.20  which  corresponds  to  an  m 
value  of  24.9  (equation  (3)).  The  cdf  values  at  1  dB  intervals  of  Xo  about 
the  median  were  also  calculated  from  the  digital  data  for  this  period  using 
equation  (2).  The  points  obtained  are  plotted  on  figure  3(b).  The  line 
on  figure  3(b)  is  obtained  from  equation  (6).  The  m-distribution 
representation  of  the  cdf  (equation  (6))  obviously  fits  the  experimental 
data  well  at  this  value  of  m. 

Figures  4  to  7  inclusive  display  similar  information  to  that  just  described 
for  figure  3,  except  the  scintillations  involved  are  characterised  by 
scintillation  indices  of  0.37  to  1.08  or  alternatively  by  m  values  of  7.2  to 
0.9  (vide  Table  1).  In  each  case  the  m-distribution  representation  of  the 

cdf  (equation  (6))  is  seen  to  give  a  good  representation  of  the  experimental 
cdf  values.  Any  discrepencies  which  occur  do  so  at  the  largest  or  smallest 
values  of  Xo  investigated,  where  the  statistics  involved  in  the  evaluation 
of  equation  (2)  are  likely  to  be  inaccurate  because  of  the  small  numbers  of 
events  encountered.  Therefore,  the  m-distribution  representation  of  the 
cdf  (equation  (6))  may  be  taken  as  a  valid  representation  of  the  actual  cdf 
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values  observed  regardless  of  the  value  of  m. 
reached  with  regards  to  the  estimation  of  bit 
(8))  given  the  value  of  S4  or  m. 


A  similar  conclusion  can  be 
error  rates  (equations  (7)  and 


4.  LONG  DURATION  SIGNALS 

The  cumulative  amplitude  probability  distribution  function  (cdf ) ,  described 
in  Section  3.1,  besides  giving  information  on  the  probability  of  reception  of 
short  period  signals  propagating  in  a  scintillating  channel,  also  describes  the 
amplitudes  of  the  fades  which  affect  signals  of  longer  period  (e.g.  >0.1  s) 
when  transmitted  via  such  a  channel .  However,  the  full  statistical  description 
of  the  effect  of  scintillations  on  long  duration  signals  also  requires  inform¬ 
ation  on  the  fading  rate.  Such  information  can  be  obtained  by  employing  either 
spectral  techniques  or  level-crossing  techniques.  These  techniques  will  be 
considered  in  the  following  Sections. 

4.1  The  autocorrelation  function 

The  autocorrelation  function  is  a  way  of  characterizing  the  rate  of  the 
scintillation  fading.  It  is  defined  in  terms  of  the  signal  intensity  I  for 
a  time  lag  r  as  follows. 


P  (T  )  =  <1(10  I(t  +  r)>/<!2  (t)> 


Here  <  >  denotes  the  average  over  the  interval  examined  (T) . 

Autocorrelograms  for  lags  (r)  up  to  50  s  have  been  calculated  for  each 
of  the  30  mm  periods  analysed  for  cdf  (figures  3  to  7).  These  appear  as 
part  (a)  of  each  of  the  figures  8  to  12  inclusive.  In  each  case  increasing 
the  lag  r  causes  the  autocorrelation  coefficient  to  fall  from  unity  at 
t  -  0,  to  zero  correlation  at  a  value  of  r  which  varies  from  correlogram  to 
correlogram.  Further  increases  of  r  sees  the  correlation  coefficient  take 
on  negative  values  and  then  return  to  oscillate  about  a  coefficient  of  zero. 
The  time  lag  at  which  the  correlation  coefficient  falls  to  0.5  (r  )  was 

measured  off  each  correlogram  and  this  measure  of  the  correlogram  width  was 
tabulated  in  Table  1. 

Perusal  of  the  values  of  tq  in  Table  1  shows  that  they  can  vary  by  a 

factor  of  three.  They  are  not  obviously  dependent  on  the  value  of 
scintillation  index  (or  m)  but  are  a  measure  of  the  average  period  of  the 
fluctuations  of  amplitude  of  the  scintillating  signals. 

4.2  Power  spectra 

The  power  spectrum  of  a  scintillating  signal  is  a  graphic  way  of  depicting 
the  frequency  structure  of  its  amplitude  fluctuations.  The  spectrum  is  the 
Fourier  transform  of  the  autocorrelation  function(ref . 17)  and  hence  is 
represented  by 


F  (f) 


1 

2tt 


27rifr 


P  (t  )  dr 


Power  spectra  for  each  of  the  half-hour  samples  listed  in  Table  1  have 
been  computed  and  these  are  displayed  as  part  (b)  of  each  of  the  figures  8  to 
12  inclusive.  The  horizontal  arm  of  each  of  the  three  small  crosses  on 
these  diagrams  indicate  the  frequency  resolution  applicable  to  that  part  of 
the  spectrum,  while  the  vertical  arm  indicates  the  magnitude  of  the  statistical 
noise  expected(ref . 17) . 
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In  general  the  spectral  shape  consists  of  a  relatively  flat  low-frequency 
spectrum  and  a  high  frequency  roll-off  of  substantially  constant  slope.  The 
break  point  between  the  low  and  the  high  frequency  slopes,  called  the  cut-off 
frequency  (f^) ,  is  a  measure  of  the  bandwidth  of  the  spectrum.  It  is  more 

clearly  defined  for  those  cases  with  a  scintillation  index  of  less  than 
0-7  (figures  8(b),  9(b)  and  10(b))  as  compared  to  those  cases  with  higher 
scintillation  index  (figures  11(b)  and  12(b)).  Estimates  of  the  cut-off 
frequency  are  tabulated  in  Table  1.  The  figures  for  the  last  two  samples 
are  enclosed  in  brackets  to  indicate  their  tentative  nature  which  results 
from  the  ill -defined  break  points  involved. 

Because  the  power  spectrum  is  the  Fourier  transform  of  the  autocorrelation 
function,  the  bandwidth  of  the  spectrum  would  be  expected  to  be  inversely 

proportional  to  the  width  of  the  autocorrelation  function.  Indeed  r  f 

0.5  c 

is  substantially  constant  for  the  first  four  samples  listed  in  Table  1. 

On  the  basis  of  this  relationship  a  cut-off  frequency  of  about  0.05  Hz  would 
be  expected  for  the  sample  commencing  at  2000  hours  on  day  236.  This  value 
lies  midway  between  the  low  and  the  high  frequency  regions  of  constant  slope 
(figure  12(b))  and  may  well  be  a  more  realistic  measure  of  the  bandwidth 
than  the  0.02  Hz  quoted  in  Table  1. 

As  illustrated  in  figures  8  to  12  inclusive,  the  high  frequency  roll-off 

in  each  of  the  spectra  involves  a  slope  proportional  to  f"P,  where  p  has  a 
value  between  2  and  5.  This  result  is  consistent  with  observations  at  other 
frequencies (ref . 18) .  It  is  obvious  from  the  value  of  p  for  each  of  the 

spectra  (Table  1)  that  p  tends  to  increase  with  scintillation  index. 

4.3  Fade  durations 

The  techniques  described  in  Sections  4.1  and  4.2,  i.e.  the  autocorrelation 
and  the  power  spectrum  techniques,  undoubtedly  give  an  insight  into  the 
temporal  aspects  of  the  full  statistical  description  of  a  scintillating 
communications  channel.  However,  they  do  little  to  alleviate  the  difficulties 
which  the  communications  circuit  designer  experiences  in  this  area. 
Consequently,  in  this  and  the  next  Section,  level-crossing  techniques  will 
be  investigated  to  ascertain  whether  they  can  offer  a  solution  to  the  problem. 

The  level-crossing  techniques  seek  to  analyse  the  scintillating  channel  by 
gathering  statistics  on  the  distribution  of  fades  at  various  signal  levels. 

In  the  first  of  these  techniques  to  be  described,  the  durations  of  fades 
below  a  particular  level  are  determined  and  tabulated  against  time  interval. 

The  levels  or  fade  margins  considered  range,  in  2  dB  steps,  from  -2  dB  to 
-10  dB  relative  to  the  median.  From  these  tabulations,  the  cumulative 
distributions  of  fade  duration  for  each  fade  margin  can  be  constructed (ref . 16) . 
In  parts  (c)  of  each  of  figures  8  to  12  inclusive,  these  cumulative 
distributions  of  fade  duration  are  plotted  for  the  signal  samples  listed  in 
Table  1.  The  points  on  these  diagrams  represent  the  results  for  each  of  the 
time  interval  cells  chosen.  The  broken  line  curves  were  drawn  freehand 
through  the  points  to  give  an  indication  of  the  likely  smoothed  distribution 
for  each  fade  margin.  The  distributions  show  that  both  the  number  of  fade 
margins  affected  and  the  number  of  fades  below  each  margin  increase  with 
scintillation  index.  This  is  the  expected  result  as  the  overall  depth  of 
fading  increase  with  scintillation  index. 

The  S-like  shape  of  the  cumulative  distribution  curves  suggests  that  the 
corresponding  simple  distribution  curves  should  be  bell  shaped.  Replotting 
the  data  as  simple  distributions  verifies  this.  However,  the  number  of 
events  per  cell  are  such  that  there  is  considerable  statistical  noise  and  no 
firm  conclusion  can  be  reached  as  to  the  nature  of  the  distributions 
(e.g.  whether  they  are  normal  distributions  or  not).  One  thing  is  clear 
however,  the  fade  duration  distributions  are  contained  between  upper  and 
lower  limits.  The  lower  limit  is  about  0.2  to  0.5  s  and  the  upper  limit 
lies  between  3  and  10  s.  There  is  some  variability  of  the  limits  within 
each  set  of  distribution  curves.  Generally  speaking  however,  the  larger 
the  fade  margin,  the  closer  the  limits.  Also  the  distribution  limits  vary 
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from  one  sample  to  another.  For  example,  both  the  upper  and  lower  limits 
of  the  distributions  for  the  sample  of  scintillations  represented  in 
figure  11(c)  are  noticeably  higher  than  those  for  the  other  samples 
illustrated.  This  is  not  unexpected  as  the  correlogram  width  for  this 
sample  is  larger  than  for  the  others  (Table  1) . 

While  this  type  of  presentation  of  the  level-crossing  data  gives  some 
insight  into  the  temporal  character  of  the  scintillation  statistics,  it, 
like  the  autocorrelation  and  spectral  approaches,  is  not  easily  related  to 
the  engineering  problem.  However,  the  next  Section  shows  how  this  data, 
expressed  in  a  slightly  different  way,  can  yield  answers  to  the  engineer's 
queries. 

4.4  Message  reliability 

Refering  to  figure  1,  it  can  be  seen  that  for  a  signal  of  a  particular 
duration  to  be  reliably  received  by  a  system  with  a  fade  margin  corresponding 
to  Xo ,  it  must  fit  completely  within  an  enhancement  of  the  signal  above 
Xo  (e.g.  within  5ts).  Consequently,  a  quantity,  message  reliability 

R(r,  Xo),  can  be  defined (ref . 16)  for  a  sample  of  scintillations  of  duration 
T  seconds,  as  the  ratio  of  the  number  of  time  intervals  or  messages  (of 
duration  r  seconds)  that  completely  fit  within  those  signal  enhancements, 
which  exceed  a  specified  signal  level  (Xo ) ,  to  the  total  number  which  could 
fit  within  T.  That  is 


R(r,Xo) 


(9) 


where  the  operator  |  |  produces  truncation  of  any  fractional  part  of  its 
argument.  Equation  (i)  can  be  rewritten  as 


T 

^  r|  5ts/r| 
R(r,Xo)  =  ~ - x 


T 


T 


Now  as  r  approaches  zero 


r|  5ts/r|  approaches 


Recalling 


o  o 

equation  (2),  it  is  clear  that  in  the  limit  r  =  0,  R(r,  Xo )  is  equivalent  to 
the  cumulative  amplitude  probability  distribution  function  P(Xo).  That  is 
the  cdf  is  the  ultimate  message  reliability  which  can  be  obtained  with  a 
scintillating  channel  for  a  particular  fade  margin  and  this  is  achieved  only 
with  very  short  duration  signals.  If  message  reliability  could  be  given 
some  analytical  expression,  it  would  provide  a  means  of  estimating  the 
increase  in  fade  margin  which  is  required  by  a  message  of  finite  duration 
over  the  value  specified  by  the  cdf  for  infinitesimal  signals,  in  order  to 
maintain  a  given  probability  of  perfect  reception. 
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4.4.1  Experimental  results 

The  second  level-crossing  technique  which  has  been  applied  to  the 
scintillation  data  involves  evaluating  message  reliability  for  message 
lengths  (r)  ranging  from  0.1  s  to  100  s  and  for  fade  margins  (Xo) 
ranging  from  -2  dB  to  -10  dB.  The  results  of  this  analysis  for  the 
data  samples  listed  in  Table  1  are  illustrated  in  part  (d)  of  each  of 
the  figures  8  to  12  inclusive.  Here  message  reliabilities  (as 
percentages)  are  plotted  against  message  length. 

In  these  figures  the  family  of  points  for  each  sample  of  scintill¬ 
ation  and  for  each  value  of  Xo  appear  to  lie  along  a  well  defined 
curve.  These  curves,  which  are  substantially  horizontal  for  the 
smaller  message  lengths,  roll-off  to  smaller  message  reliabilities  at 
the  longer  message  lengths.  A  check  shows  that  for  each  of  the  curves, 
the  message  reliability  found  for  the  smallest  message  length  examined 
(0.1  s)  corresponds  to  the  cdf  value  expected  for  the  values  of  m  and 
Xo  involved  (figure  2) .  Indeed,  the  departure  from  the  cdf  value 
does  not  become  appreciable  until  the  message  length  exceeds  about 
5  s. 

Increasing  the  scintillation  index  results  in  both  a  general 
decrease  in  message  reliability  for  any  fade  margin  and  also  increases 
the  number  of  fade  margins  so  affected.  Also  comparison  of 
figure  11(d)  with  figures  10(d)  and  12(d)  shows  that,  while  they  each 
represent  severe  scintillation,  the  high  message  length  roll-off  of 
message  reliability  occurs  at  a  larger  value  of  r  in  figure  11(d)  than 
in  the  other  two  figures.  This  is  consistent  with  the  larger  auto- 
correlogram  width  (r ^  ,.)  for  the  data  sample  from  which  figure  11(d) 

is  produced  compared  to  the  values  of  tq  $  for  the  data  samples 

corresponding  to  the  other  two  figures  (Table  1).  These  are,  of 
course,  results  which  are  expected  on  the  basis  of  the  discussions  of 
Sections  4.1  to  4.3  inclusive.  They  reaffirm  the  suggestion  made 
earlier  that  a  statistic,  which  purports  to  fully  represent  the  effects 
of  scintillation  on  a  communications  channel,  must  take  both  the 
amplitude  and  the  time  dependencies  into  account. 

4.4.2  Empirical  theory 

If  message  reliability  is  to  be  developed  as  a  quantity  which  is 
useful  in  the  engineering  of  scintillating  communications  channels, 
it  is  necessary  to  give  it  analytic  expression.  This  Section  will 
develop  such  an  expression  in  a  semi-empirical  way. 

The  probability  of  the  disturbed  signal  being  at  a  level  X  (pr  (X)) 
is  known  (equation  (4)).  The  message  reliability  depends  on  this 
probability  and  also  on  the  probability  of  the  disturbed  signal  at  or 
above  a  level  Xo  remaining  at  or  above  that  level  for  a  time  r. 

If  this  latter  probability  is  called  the  abiding  probability 
a(Xo,r),  then  the  message  reliability  can  be  written 


R  (Xo  >  t0  )  = 


p'  (X)  a(Xo,r)  drdX 


Xo  J  r0 


P  (X)  dX  f  a(Xo,r)  dr 
Xo  r0 
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Using  equation  (6)  this  becomes 


R(Xo  ,r0) 


P(Xo) 


a (Xo ,r)  dr 


r0 


or 


R(Xo  ,r0 )  =  P  (Xo )  A(Xo,ro)  (10) 

That  is,  the  message  reliability  is  equal  to  the  cdf  times  a 
reliability  factor  A(Xojr0)  where 

T 

A  (Xo  ,  r0 )  =  I  a(Xo,r)  dr 

•>  r0 


It  has  been  demonstrated  experimentally  (Section  4.4.1)  that 
A(Xo>r0)  is  equal  to  or  less  than  unity.  It  approaches  unity  as  r0 
approaches  zero  and  its  departure  from  unity  is  only  appreciable  for 
values  of  ro  in  excess  of  about  5  s. 

If  the  data  of  each  part  (d)  of  figures  8  to  12  inclusive  are  re¬ 
plotted  with  linear  message  length  scales,  the  groups  of  experimental 
points  fall  along  straight  lines  of  negative  slope  which  intercept  the 
message  reliability  axis  at  the  appropriate  cdf  values.  That  is 

logio  R(Xo,ro)  =  logio  P(Xo)  -  kr0 


or 


R(Xo,r0)  =  P(x0)  10  kro  (11) 

It  is  then  found  that  a  log- linear  plot  of  Xo  versus  k  yields  a 
straight  line  of  negative  slope  for  each  of  the  data  samples  represented 
in  figures  8  to  12  inclusive.  Hence 

logio  Xo  =  logio  L  -  ck 

where  L  is  the  intercept  on  the  Xo  axis.  This  can  be  rewritten  as 

10k  =  (L/Xo)1/C 

which  when  combined  with  equation  (11)  yields 

R(Xo,r0)  =  P  (Xo )  (L/Xo)"r°/c  (12) 


Since  the  quantity  L  is  related  to  signal  level,  it  is  most  likely 
to  be  dependent  on  scintillation  index.  Indeed  it  is  found  that,  for 
the  five  scintillation  samples  listed  in  Table  1,  L  is  directly 
proportional  to  scintillation  index,  the  relationship  being 
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L  =  -17.3  S4 


(13) 


Since  2M  =  17.37  (equation  (5))  and  m  =  1/S4  (equation  (3)), 
equation  (13)  can  be  conveniently  restated  as 


L  =  -2M/Vm 


(14) 


In  a  similar  way,  since  c  is  related  to  message  length  ro ,  c  is  most 
likely  to  be  related  to  the  time  dependence  of  the  scintillation.  It 
is  found  that  for  the  five  samples  of  scintillation  illustrated,  c  is 
directly  proportional  to  r  the  relationship  being 


c 


10.4  t 


0.5 


(15) 


Combination  of  equations  (12),  (14)  and  (15)  gives  finally 


-0.096  t0/t 

R(Xo,r0)  =  P(Xo)  (-2M/vfiTxo) 


or 


R(Xo  ,r0) 


(-2M/\4n  Xo) 


u .  5 


M  T(m) 


r0 


exp[  m  [  2x/M-exp  (2x/m)j  ]  dx 

Comparing  equation  (16)  with  equation  (10),  it  is  seen  that 


(16) 


(17) 


A(Xo  ,r0) 


a(Xo,r)  dr 


ro 


-0.096  r0/r 

(-2M/yfiT  Xo  ) 


(18) 


Consequently,  the  abiding  probability  a(Xo,r)  is  given  by 


a(Xo  ,r) 


0.096 


0.5 


-0.096  r/T 

l°g(-2M/yfii  Xo)  ( - 2M/y^n  Xo) 


To  further  demonstrate  the  validity  of  equation  (17),  it  was  used  to 
simulate  the  variation  of  message  reliability  with  message  length  for 
the  conditions  pertaining  to  each  of  the  families  of  experimental 
points  plotted  on  parts  (d)  of  each  of  the  figures  8  to  12  inclusive. 
The  broken  line  curves  on  these  figures  are  the  result  of  these 
simulations  and,  as  can  be  seen,  good  agreement  is  obtained  with  the 
experimental  points. 
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4.4.3  Circuit  design  chart 

The  above  discussion  indicates  that  message  reliability  is  a 
convenient  single  parameter  for  characterizing  the  effects  of 
scintillations  on  a  transionospheric  communications  channel.  For 
signals  of  short  duration,  message  reliability  is  equivalent  to  the  cdf 
which  is  a  measure  of  circuit  performance  already  in  common  use 
(ref . 11, 16, 18) .  An  algorithm  for  the  computation  of  the  cdf  can  be 
based  on  equation  (6) ,  while  for  most  communications  circuit  design 
purposes,  the  graph  presented  as  figure  2  would  provide  cdf  values  to 
sufficient  accuracy. 

For  signals  of  duration  greater  than  about  a  second,  message 
reliability  takes  on  values  less  than  the  cdf  value,  the  reliability 
factor  A(Xo,r0)  (equation  (18))  being  the  appropriate  multiplying 
factor.  From  equation  (16)  this  is 


-0.096  ro /T  0 _  5 

A  (Xo  >  ro )  =  R(Xo,r0)/P(Xo)  =  (-2M/vfiTxo)  '  (19) 

Equation  (19)  could  be  used  as  the  basis  for  the  construction  of  an 
algorithm  for  the  calculation  of  the  reliability  factor.  An  algorithm 
for  the  complete  computation  of  message  reliability  (R(Xo,r0))  would  be 
based  on  both  equations  (6)  and  (19) . 

Examination  of  the  form  of  equation  (19)  suggests  that  a  plot  of 
A(Xo,r0)  versus  r/r  and  parametric  in  XoV^n  would  provide  a  useful 

chart  for  circuit  design  purposes.  Figure  13  is  such  a  diagram. 

Note 


A  (Xo ,  ro )  =  1  for  Xo  V®"  =  -2M  =  -17.37 

A  (Xo  >  ro )  <  1  for  Xo  yfiT  <  -2M 


while  values  of  Xo  and  m  which  make  Xo  >  -2M  are  physically 

incompatable.  . 

Figures  2  and  13  allow  message  reliability  to  be  determined  quickly 
for  particular  fade  margins  Xo  and  message  lengths  ro  given  a  knowledge 
of  m  (or  scintillation  index)  and  rQ  5* 


5.  DISCUSSION  AND  CONCLUSIONS 

In  order  to  facilitate  the  design  and  to  evaluate  the  performance  of  satellite 

communication  links  which  are  subject  to  scintillation,  it  is  necessary  to 

characterize  scintillation  in  a  manner  different  to  that  which  has  been  found 
valuable  when  using  scintillation  for  the  study  of  ionospheric  irregularities. 

The  engineer  designing  a  satellite  communications  link,  needs  to  know  what 
fraction  of  the  time  communications  are  likely  to  be  disrupted  and  how,  by 
suitably  choosing  certain  parameters  of  the  system,  he  can  minimise  this  loss. 

The  cumulative  amplitude  probability  distribution  function  (cdf)  has  been  used 
as  such  a  design  parameter.  However,  it  has  been  shown  that  this  is  appropriate 
only  for  signals  of  short  duration  (<  1  s)  and  then  strictly  only  for  amplitude 
modulated  signals. 

If  signals  with  durations  of  a  second  or  more  are  to  be  catered  for,  the  cdt 

is  no  longer  a  suitable  design  parameter  and  account  has  to  be  taken  of  the 

temporal  character  of  the  scintillations.  The  quantity  message  reliability 
which  is  defined  above,  provides  to  be  a  parameter  which  is  suitable  m  this 
regard.  It  is  valid  for  signals  of  all  durations  and  includes  the  cdf  estimate 
as  a  special  case  when  the  signal  duration  is  small  (<  1  s) .  Further,  message 
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reliability  has  been  shown  to  be  a  meaningful  statistic  for  all  levels  of 
sctnti nation  activity.  However,  like  the  cdf,  it  is  strictly  valid  only  for 

amplitude  modulated  signals,  though  the  approximation  involved  for  other  types 
of  modulation  improves  with  increasing  message  duration. 

The  relationship  between  message  reliability  and  cdf  has  been  investigated 
in  a  semi -empirical  way  using  UHF  scintillation  data.  The  formula  so  developed 

aiiows  message  reliability  to  be  evaluated  for  any  fade  margin  (Xo)  and  message 


0.5 


) 


length  (r)  given  the  m  factor  and  the  half  width  of  the  autocorrelogram  ( t 

appropriate  to  the  scintillating  channel.  A  design  chart  has  been  produced 
which  gives  expression  to  this  formula  and  allows  message  reliability  to  be 
evaluated  readily  in  terms  of  r/rQ  g  and  Xov^- 

The  data  used  in  the  analysis  leading  to  the  formulation  of  the  message 
reliability  expression  applies  to  one  carrier  frequency  and  one  location.  The 
formula's  validity  for  other  frequencies  and  locations  obviously  need  to  be 
evaluated.  Also  bit  error  rate  measurements  obtained  from  phase  modulated 
signals  are  needed  to  further  evaluate  the  usefulness  of  the  cdf  and  message 
reliability  concepts  for  such  signals. 
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(b)  cumulative  amplitude  probability  distribution 


Figure  3.  Thirty  minutes  of  scintillation  (S4  =  0.20) 
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Figure  4.  Thirty  minutes  of  scintillation  (S4  =  0.371 
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Figure  5.  Thirty  minutes  of  scintillation  (S4  =  0.65) 
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(b)  cumulative  amplitude  probability  distribution 


Figure  6.  Thirty  minutes  of  scintillation  (,S4  =  0.86) 


ERL-0074-TR 
Figure  7 


(a)  analog  record 


oCUMULFITIVE  PRQB.  DIST.  DAY236  20  0  S4rl.08  Mr  0.9g 
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Figure  7.  Thirty  minutes  of  scintillation  (S4  =  1.08) 
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